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Abstract: Control of spatial quantum correlations in bi-photons is one of 
the fundamental principles of Quantum Imaging. Up to now, experiments 
have been restricted to controlling the state of a single bi-photon, by 
using linear optical elements. In this work we demonstrate experimental 
control of quantum correlations in a four-photon state comprised of two 
pairs of photons. Our scheme is based on a high-efficiency parametric 
down-conversion source coupled to a double slit by a variable linear optical 
setup, in order to obtain spatially encoded qubits. Both entangled and 
separable pairs have been obtained, by altering experimental parameters. 
We show how the correlations influence both the interference and diffraction 
on the double slit. 
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1. Introduction 

Works exploring the quantum effects on multi-photon spatial interference and images have 
traditionally focused on experiments using a pair of photons IBQ. Other important results 
obtained in the past decades, involving the spatial degree of freedom of photon pairs, were the 
complementarity between entanglement and single-photon coherence |]6l, increased resolution 
on imaging when using entangled photons 17t-(9j|, the spatial quantum eraser ifTOl . measurement 
of the two-photon de Broglie wavelength 11 1II12I . and the experiments encoding logical qubits 
and qudits onto the path state of photons going through slits ITsm sl . 

In this paper we demonstrate control over the quantum correlations of two biphotons that are 
encoded in a four qubit space by using a double slit. Our setup is able to generate two pairs of 
highly entangled photons in one configuration and four nearly-separable photons in another We 
further present a way of characterizing the state generated, without having to resort to fourfold 
coincidence detection, by exploiting the fact that each pair contains one photon of horizontal 
polarization and one of vertical polarization. This allows us to observe both the contribution 



that is common to single and double pair spontaneous parametric down-conversion (SPDC) 
events and the one that is exclusive to double pair events, i.e. coincidences between photons 
with the same polarization. Finally, we also observed that the diffraction behavior depended 
on the quantum correlations of the photons, producing an asymmetrical pattern in the case 
of entangled photons and the expected diffraction pattern for the case of photons in a nearly- 
separable spatial state. 

By encoding the transverse profile of photons on a discrete space with the use of slits, it 
becomes simple to scale the dimensionality, from qubits to qudits, by using more slits. This 
experiment, then, opens a path for future works using multi-pair qudit states. As it is custom- 
ary, double-slits and multi-slits are characterized by their width 2a and by the center-to-center 
distance 2d. 

2. The four-photon state 

In order to exploit the possibility of scaling experiments to include quantum systems of dimen- 
sions greater than 2, it has become usual to encode logical qubits and qudits on the transverse 
momentum and/or position state of photons. Our experiment encodes the logical qubits |0) and 
1 1) on the condition that the photon passes through the left or right slit, respectively. 

As presented for the two-photon case in 01611181 . the most general state of a pair of qubits 
generated using a double slit and assuming a symmetrical illumination can be written as: 



where and |0+) are maximally entangled states of two photons A and B known as Bell 
states: 



Up to now, experiments have been restricted to controlling the spatial state of a single bi- 
photon, by using linear optical elements. In this work we demonstrate experimental control 
of quantum correlations in a four-photon state comprised of two pairs of photons. By using 
a pulsed pump and a high-efficiency periodically-poled KTiOP04 crystal |fT9]j2Tl l. the second 
order term in x^^^ on the expansion of the time-evolution operator associated to the SPDC 
becomes relevant. In our experimental conditions, this term results in the generation of two 
distinguishable pairs of photons. 

The pump is a 413 nm pulsed beam obtained by second harmonic generation on a BBO crys- 
tal pumped by an 826 nm femtosecond pulsed laser that has pulses of 200 fs and a repetition 
rate of 76 MHz. The mean power of the pump beam when falling on the crystal is of the order 
of 80 mW. It is loosely focused by a 30 cm spherical lens onto a 10 mm long, type-II PPKTP 
grown for collinear and degenerate SPDC centered on A = 826 nm. The created pair will con- 
sist of one horizontally (H) and one vertically (V) polarized photon. Polarization was used as 
an auxiliary mode, in order to separate the photons of each pair into two distinct physical sub- 
systems A and B, much in the same way non-collinear experiments that encode the qubits using 
polarization use the down-conversion linear momentum modes. 

The general state, considering one or two pairs cross the double-slit is, up to appropriate 
normalization Il22l . 



\^>H,v) = cos(a/2) \w^)+ e'f sin(a/2) 
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where 77 is the generation efficiency, depending on both the non-linear susceptibihty and the 
pump intensity, and M is a normalization constant. 
The four-photon component can then be written as 

|1'//.y)i®|^//,y)2= (cos(a/2)|v/+)+e"Psin(a/2)|(/)+))®' (5) 

Since there is no genuine multipartite entanglement, it is sufficient to characterize the entan- 
glement by the concurrence ||23l of the single pair, which is given by II16II18I 

= y^l-(cos(^)sin(a))2 (6) 

States with ^ = ;r/2 or a = 0, tt are maximally entangled. Conversely, states with ^ = and 
a = ;r/2 are separable. 

The presence of two photons of horizontal and vertical polarization also gives us the opportu- 
nity to consider the partial state of same-polarization photons. In practical terms, this is obtained 
by using a polarized beamsplitter to create two detection branches, one for horizontally and one 
for vertically polarized photons. 

Since the two pairs are not entangled, the purity of the partial state, defined as T^J^iPpaniai) 
and given below, will be used as a relevant indicator on those same-polarization coincidences, 
as it is directly related to the degree of entanglement in each pair: 

i32> = i(l + (cos((p)sin(a))2)l (7) 

The combination of Eq. ^ and Eq. (|7| allows us to write 

^ = ^{2-'^y. (8) 

3. Control and characterization of spatial quantum correlations 

In order to control the quantum correlations in the four photons, we will resort to a tunable 
linear optical scheme. Here we present two extremal conditions, that is, the pairs of photons 
coming out of the double slit are either highly entangled or separable. The correlations after 
the double-slit depend on how the bi-photon field generated inside the crystal is projected onto 
the plane of the apertures. Influence over the state encoded on the slits can be exerted either 
by controlling the pump field, the propagation of the bi-photons from the crystal to the slits 
or both. Since the bi-photon field is not affected in one of the transverse dimensions by the 
double-slit (the larger dimension, i.e. the y direction), we restrict our analysis to the x direction 
(the smallest slit dimension). 

More specifically, the two-photon field in the crystal can be written, in the transverse momen- 
tum variables, as III8II24I 

^iqi,q2)^Ep{qi+q2)^iqi-q2), (9) 

where {^1,^2} are the transverse components of the momentum of photon 1 and 2 in the bi- 
photon. The tilde indicates that the functions refer to the transverse momentum space. Therefore 
Ep is the angular spectrum of the pump beam and ^ is the phase-matching condition for the 
PPKTP. 

In most cases, including our experiment, Ep is the angular spectrum of a gaussian beam. 
Hence, Ep{qi +^2) °^ -^(6xp( — (^It^)^)) when the beam waist, Wp, is in the center of the 
crystal, with ^ corresponding to the Fourier transform and xi and to the transverse position 



of photons 1 and 2. If the last condition is not satisfied, phase components due to the curvature 
of the wavefront must be considered as well. 

For our experimental conditions, the phase matching function is 



^ cx: sine ( ^0 



iHeffCOsPDc/c) 



(10) 



where 0o is a phase-mismatch term that we consider to be negligible, COspdc is the frequency of 
the degenerate down-converted photons and the effective refractive index itgff = 2(«i«2)/(«i + 
ni) takes into account the fact that the two photons in type-II SPDC propagate differently inside 
a crystal of length L. A more thorough description of the phase-matching function can be found 
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Fig. 1 Experimental setup for control of transverse correlations over the double slit (DS). Cou- 
pling is done so that either the position (a) or momentum (b) degree of freedom of the bi-photons 
is mapped onto the double slit plane |26| . The first uses a cyhndrical lens (CL) of fi = 5 cm and 
a spherical lens (SL2) of f2 = 20 cm to project a magnified image of the crystal center, while 
the latter uses only the = 20 cm lens to project the Fourier transfonn of the same plane onto 
the plane of the double slit. The slits are 80 |jm wide and have a separation 2d = 240 \m\. SLl 
is the / = 3 1 cm spherical lens that focuses the pump beam on the PPKTP crystal and DM 
is a dichroic mirror that reflects the 413 nm beam while allowing the 826 nm down-converted 
photons to pass. 



Figure[T]shows the experimental setups used in this work to control the quantum correlations 
of the bi-photons at the double-sUt plane. We adopted two configurations: |l(a)| to generate 
highly entangled states and |l(b)| to prepare separable states. 

Scheme [l (a)| projects a magnified image of the crystal plane onto the double-slit. This ge- 
ometry is achieved by adopting two lenses, a cyhnder lens with focal lens /i = 5 cm and a 
spherical one with /2 = 20 cm, so the magnification factor of the image reads M = /2//1 =4. 



We adopt a cylinder lens in order to create a crystal image in the direction perpendicular to the 
slits, the x-direction, while keeping the far-field image in the y orthogonal direction. By such 
approach the dominant contribution to the bi-photon amplitude transmitted by the double-slit is 
the one associated to the phase-matching function. That is so because, in the magnified image 
of the crystal center projected on the slits, the gaussian profile of the pump component is much 
wider than the Fourier transform of the sine function of the phase-matching condition and can 
be considered constant. 

Configuration 1 1 (b)| relies on the far-field imaging onto the double slit. This scheme is ob- 
tained by adopting a single spherical lens with focal length / = 20 cm placed at a distance / 
from both the center of the crystal and the plane of the slits. With this geometry the bi-photon 
amplitude is almost constant along the double-slit apertures leading to the generation of separa- 
ble states. 

A detailed discussion of the relation between the parameters a and <p of Eq. ([U and the 
adopted geometries can be found in ITSl . in the form of what was called an engineering param- 
eter p: 

p = exp(/(p)tan(a/2) = -^^^, (11) 

that gives a clear method for estimating the values of a and (p by knowing the form of the two- 
photon amplitude function A(xi ,X2) specifically when xi = X2 ~ d and xi = —X2 = d, where 
xi and X2 are coordinates on the plane of the double-slit and d, as previously mentioned, is 
the distance from the center of one of the slits to the center of the slit pattern. In this case, the 
distance between the center of the slits is 2d. We refer to the values calculated through that 
method when giving expected values for a and (p. For our experimental setups, the expected 
values are a„ = 173° and (p„ = 180° in the configuration of Fig. |l(a)| By changing to the 
configuration of Fig. [T(b)| we expect at, ~ 85° and (ph = 0°. 

Let us now discuss the measurements necessary to characterize the states generated by the 
setups presented in Fig.[T] 

The parameter a can be inferred by measurements on the image of the slits, which corre- 
sponds to the {|0), |1)} basis. Measurements for determining (p, however, require us to go to 
the diagonal basis. For qubits encoded on transverse spatial variables, this can be achieved by 
looking at coincidence measurements on the pattern generated in the far field of the slits, that 
gives information on both a and (p. One way to access the far field is to use a lens with focal 
length ffi in the /-/ configuration. This maps the Fourier transform of the slits onto positions 
on the detection plane, as shown in Eq. 5-19 of [26 |. 

The probability of coincident detection of two photons from the same pair, one of horizontal 
polarization and the other of vertical polarization, when using the far-field strategy, can be 
derived as a function of the position of both detectors involved, P{xi,X2). To calculate it, we 
first assume the contribution of only a single pair of photons and then follow Mandel and 
Wolf EH to obtain 

PixuX2)(H.V)'^\E+ixi)E+{x2)\1'H.v)\\ (12) 

where E^ixi), for j = {H,V) and / = (1,2) is calculated by propagating the electrical field 
operators from the plane of the slits to the detection plane, taking into account the lens placed 
at a distance fpi from both planes. This, as mentioned before, means taking the optical Fourier 
transform of the product of the aperture function and the amplitude of the bi-photon on the 
plane of the slits ifTSlfTSl . 



After performing the calculations, we obtain the explicit form of P{xi,X2)(^H,vy 



P{xi,x2)(H.v) °^ 1 +cos^(a/2)cos()3(xi -X2))+ 

I sin(a) I cos((j!)) [cos(j3 (xi )) + cos()3 (x2))] + 
sin^(a/2)cos(j3(jci +X2)). 



(13) 



This function can be represented on a plane, which we will call the coincidence map. The 
same approach can be adopted to represent the coincidence counts detected in the experiment, 
allowing us to compare experimental and theoretical maps. We see that it depends of both a 
and (p, meaning the coincidence maps can be used to characterize both parameters of the state. 

The following step is to characterize the measurements sensitive exclusively to the two-pair 
contribution \^Hy)j'^ \^h,v)ji- 

Since we work with two pairs of photons, the most natural type of measurement would be 
fourfold coincidence detection. This, however, is impractical due to the fact that only a small 
fraction of the total signal is collected by the detectors when in the far-field regime. As four- 
fold coincidences depend on the overall efficiency of four detectors, it can require extremely 
long detection times. To circumvent this predicament, we use the fact that our source of pho- 
ton pairs is a type-II crystal and, thus, coincidences on the same polarization branch can only 
come from down-conversion events in which two pairs were created. Another factor that favors 
same-branch coincidences as a method for investigating the two-pair signature is that only one 
beamsplitter (BS) is involved, whereas two are necessary on fourfold detection events, as can 
be seen in Fig.|2] 

These coincidences can be mapped in much the same way as those from two photons with 
orthogonal polarization. Even though double-pair events are not the major contribution, the 
probability of detecting H-H or V- V coincidences is still much larger than of detecting fourfold 
coincidences, since only two detectors are involved in the process and, therefore, only two 
detection efficiencies are taken into account. 

The expression for coincidence detection of photons on the same polarization branch is: 



in which j can be either H or V. Indexes / and // correspond to the two pairs of photons. 
The density matrix pj is given by p'j = ^ri{\^'Hv) {^hv\), with Tr,- being the partial trace over 
the photon of polarization ; 7^ j. Note that the same derivation can be used to calculate the 
probability of detecting coincidences between photons of different polarization and different 
pairs. 

Equations ( fT3] l and (fT4l i have been derived under the assumption that the incident two-photon 
field can be considered constant over the width of a single slit, which is taken to be very narrow. 
This allowed us to factor out the diffraction term as an enveloping function that is radially 
symmetric around the origin in the plane formed by xi and X2- As we will discuss in Section |5] 
this is not true for all situations, but is a valid approximation if one is interested only in studying, 
for instance, the interference behavior close to the center of the detection plane. 



P{xi,X2){(H.H);(v.v)} °= 1 + (sin(a)cos(9))2cos(xi)cos(x2) + 

sin(a)cos(9)(cos(xi) +cos(x2)). 



(14) 



To calculate it, we use 






Fig. 2 A schematic view of the detection apparatus, showing how the incoming photons are di- 
vided into two branches by a PBS. On each branch the photons are further divided by a BS. This 
allows us to make coincidence detection between two detectors on the same branch or between 
detectors on different branches. Detections on the same branch originate only from double-pair 
events, while the raw different branch detections are a sum of contributions from single-pair and 
double-pair events. In this figure, the focal plane of FL corresponds to the detection plane, upon 
which all four detectors are free to move parallel to the plane of the table and perpendicular to 
the direction of beam propagation after passing through the beamsplitters. DS is the double-slit. 



We will now discuss in detail the experimental procedure used to obtain the coincidence 
maps. The detection setup, shown in Fig.[2]is composed of a spherical lens FL of focus length 
/ = 30 cm and four avalanche photo-diodes (APD) detectors divided into two polarization 
branches, as indicated above by the dashed bounding boxes. Each branch corresponds to one of 
the outputs of a polarized beamsplitter (PBS) and contains a balanced beamsplitter (BSl and 
BS2). Thus, coincidences can be detected between APDs on the same polarization branch, to 
which we refer as SBC (Same Branch Coincidences), and on different branches, referred to 
as DBC (Different Branch Coincidences). Each detector is mounted on a manual translation 
stage in the x direction, defined as being parallel to the optical table plane and orthogonal to 
the propagation direction of the beam. Additionally, a 200 ^.m pinhole and an interference filter 
(AA = 40 nm) centered on 826 nm are placed in front of the APD. 

Coincidences are detected in a 1 ns window using a custom electronic circuit. Such a time 
frame is large enough to collect all events originated from a single 200 fs pulse but small enough 
so that there can be no coincidence between photons from two different pulses, allowing us to 
consider each pulse as an isolated SPDC event that may generate either zero, one or two pairs of 
photons. The maps are the result of 21 scans (17 in the case of Fig.|5]l where two detectors, one 
in each branch remain in a fixed position and the other two are moved in steps of 100 i_im (for 
fine scans over the main interference region) or 500 i.im for broad scans, meant for observing 
the diffraction pattern over a larger area. 

In addition to generating the maps, the data collected was also used to evaluate efficiency of 
our source. To calculate the contribution of single and double pair events, we have estimated 
the detection efficiency of our apparatus starting from the measured count and coincidence 



rates of detectors. Considering uncorrelated spatial photon pairs generated with an average 
pair number per pulse p, detected with an overall efficiency J, for an experiment with time 
duration T and repetition rate R, the average counts of the detector are equal to A^i = jpRx, 
while the measured coincidences are Nu ~ 'f'pR'Z. The value of 7 can then be obtained as 
Nn/N\. Spatial correlations between the two photons can modify the previous relation, so the 
data used corresponds to the maxima of detection on each map. 

4. Experimental results 

Following the experimental procedures detailed in the previous section, we obtain coincidence 
maps on the far field of the double slits, which will serve as fingerprints for identifying the states 
generated. The parameters obtained from these maps were then compared to those predicted by 
the theoretical model of our experiment, as presented in section[3] 

Maps on Fig. [3] El and |6] were generated from a matrix of 21x21 experimental points. Those 
on Fig. |5] originate from a 17x17 grid, although the separation between points are the same as 
that of Fig. |6l meaning that a smaller field of view was imaged, from -4 to 4 mm instead of 
-5 to 5 mm. Data points have been normalized relative to the maximal values of each map, for 
ease of comparison to the simulated maps based on the theoretical model from Eq. ( fT3] l and 
Eq. (fT4] i. Additionally, the latter are generated over the same grids of points, corresponding to 
the same values of xi and X2- This means that artifacts due to interpolation are not exclusive to 
the experimental maps and should present themselves in the same manner on the simulations, 
minimizing discrepancies that may influence the interpretation of the processed data. There 
are, however, some artifacts present only on experimental maps, due to noise on the data set. 
On the case of maps on Fig. [3] and Fig. |4] the dominant aspect is noise, since the structures 
are much larger than the spacing on the grid. On Fig. |5]and Fig. |6l however, the size of the 
finer structures is comparable to that of the spacing between points on the grid. This makes 
the aliasing contribution more significant. However, in those two figures the most interesting 
structures, those due to diffraction, are larger than the grid spacing. 

Equation ( fT3b was used to simulate the different branch coincidence map for comparing with 
the experimental DBC map. Since this equation doesn't take in account coincidences between 
photons of orthogonal polarization but from different pairs, we have subtracted that contribution 
from N12 by estimating the accidental coincidences as Nacc — Ni *N2/ (^t). These accidental 
coincidences, calculated for each measured point, were also subtracted from the maps. 

Using the corrected data, we also characterize our source. Calculating the average pair num- 
ber p ~N\/ {yRz) we obtain pa = 0.49 and pi, = 0.96, respectively for the setup of Fig. [T(a)] and 
|l(b)| We attribute the difference between the two values to the different geometries and to the 
spatial correlation observed with setup |l(a)| The effective single pass parametric gain, calcu- 
lated after the double slit, is 0.65 for setup |l(a)] and 0.87 for |l(b)| The probabilities P{n) of gen- 
erating n pairs are, for[T(a)]and[T(b)] = 0.29, ^(l)^ = 0.37, P(2)„ = 0.07, P{2)i, = 0.19, 
p{3)a = 0.01 1 and P(3),, = 0.066. 

We now examine in detail the results of the entangled state obtained by using the setup of 
Fig. |l(a)| and, after it, the near- separable state generated by the setup indicated in |l(b)| 

4.1. Entangled state 

In Fig. [3] the SBC map, originated only by multi-pair events, confirms that the partial state 
is indeed mixed. The fact that a small degree of structure, still present on the simulated SBC 
map, is not observed experimentally is attributed to detection fluctuations that mask the low 
visibility pattern. The interference pattern observed on experimental DBC maps, which reflect 
mainly the single-pair contribution, is consistent with near maximal entanglement on the qubit 
picture. Those two results are interpreted as a sign that two distinct pairs of photonic spatial 



qubits with high purity and path entanglement have been generated. 




Fig. 3 Normalized two-plioton measured and simulated coincidence maps for a four-photon path 
state presenting a high degree of spatial correlation in each bi-photon. The upper row, from (a) 
to (c), includes the measured coincidences on the H branch (D2D4). the calculated uncorrelated 
pattern corresponding to the product of the single counts (D2*D4), weighted by the repetition 
rate of the pump, and the SBC simulation, generated with Eq. il4\ . The lower row, from (d) to (f), 
corresponds respectively to the coincidences between the H and V branches (D2D3) and to the 
same map corrected by subtracting the uncorrelated contribution (D2D3-D2*D3), as well as a 
simulated DBC map, generated with Eq. The simulation parameters used are = 176° 
and 9"'' = 170°, corresponding to V 0.997 ±0.003. For a measurement time of 25 s per 
point, the maximum in map (a) coiTesponds to 549 coincidence counts and in map (d) to 1630 
coincidence counts. The measured maps were obtained at the Fourier plane of Fig.|2] after the 
photon pairs cross the double-slit. 

In Fig.[3td) and Fig.[3te) one can notice that the visibiHty increases substantially after sub- 
tracting from the total detections the coincidences between photons of different pairs generated 
by the same pulse, and transmitted through the double-slit. Our result then demonstrates that 
highly entangled pairs can be generated even in the regime where many bi-photons are produced 
by a single pump pulse. We also note that, on each of the three fringes visible in the pattern, the 
maxima are situated on xi =X2 = ±500 |xm or ixm. This period of 500 ^m coincides with the 
step used on the measurements of Fig. |5] 

As previously stated, the expected values for the two characteristic parameters for the state 
were tta = 173° and % = 180°. In order to maximize the agreement between the experimental 
data and the simulation, those two parameters were changed to Ua"'' = 176° and ^"'^ = 170°. 
Such a change for (p would be expected if the position of the pump waist and of the focal 
plane of the cyllindrical lens differed by approximately 200 microns, which is well within the 
experimental uncertainties of our setup. For the values used in the simulation, the concurrence 
is "if « 0.997 ± 0.003 and the purity for the partial state is ^ « 0.254 ± 0.003, compai-ed to a 
minimum of 0.25 for a maximally mixed state. 



4.2. Near-separable state 

Let us now analyze the case of the two-photon coincidence maps generated by the setup of Fig. 
|l(b)[ for which the theoretical prediction is a four-photon state of negligible entanglement in the 
path variables, after the double-slit. In this case, the single- and double-pair contributions can 
be considered as a tensor product of two and four distinguishable single-photon pure path states, 
respectively. This is evidenced by the fact that the experimental maps for SBC and DBC are 
similar For a state of two pure and separable pairs, the partial state detected will also be pure. 
Conversely, for maximally entangled path states the partial state is maximally mixed. Thus, we 
see no interference on the map of Fig. Ha), but the maps of Fig. Ua) and Fig. Ud) are equal. 
This is in direct relation to what was shown in Il6ll28l. 

For this state, the expected values of the two characteristic parameters were Ub ~ 85° and 
(ph = 0° . As done during our analysis of the entangled case, those two parameters were changed 
to a^*'' = 86° and (j?^^'' = 5° in the simulated maps. Such a change for both a and (p is, again, 
well within the experimental uncertainties of our setup. For the values used in the simulation, 
the concurrence is « 0.11 ±0.03, which results in a high purity for the partial state, ^ sa 
0.988 ±0.004. 
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Fig. 4 Normalized two-plioton measured and simulated coincidence maps for a four-plioton path 
state of negligible spatial correlation. The upper row, from (a) to (c), includes the measured 
coincidences on the H branch (D2D4), the calculated uncoiTelated pattern corresponding to the 
product of the single counts (D2*D4), weighted by the repetition rate of the pump, and the 
SBC simulation, generated with Eq. U4t . The lower row, from (d) to (f), corresponds respec- 
tively to the coincidences between the H and V branches (D2D3) and to the same map corrected 
by subtracting the uncorrelated contribution (D2D3-D2*D3), as well as a simulated DBC map, 
generated with Eq. l ll3t . The simulation parameters used are a"'' = 86° and ip"'' = 5°, corre- 
sponding to 'i^' « 0.1 1 ± 0.03. For a measurement time of 25 s per point, the maximum in map 
(a) corresponds to 122 counts and in map (d) to 247 counts. The measured maps were obtained 
at the Fourier plane of Fig.|2] after the photon pairs cross the double-slit. 

One should note that the fluctuations observed in Fig.|4]are due to the lower statistics, com- 
pared to those of Fig. [3] In the former, Poissonian fluctuations are more significant. Another 
important point to be observed about the maps in Fig. |4]is that the peaks are 1 mm away from 



each other, twice the step used on Fig.|6] 
5. Diffraction 

In order to fully understand our experiment, we now turn our attention briefly to the subject of 
the diffraction envelope. We present a qualitative discussion of its characteristics for the two 
extremal cases of transverse spatial quantum correlations that were generated in this work. As 
in the previous section, we find it appropriate to discuss the case of a highly entangled state and 
of a near-separable state separately. 

5. 7. Entangled state 

Figure ID shows a larger area for the same experimental conditions as Fig. [3] The asymmetric 
nature of the diffraction envelope for different branch coincidences is clear, with a larger diffrac- 
tion on the {x\ = — X2) = x_ direction and a narrower one for {x\ = X2) ~ x+, if compared to 
the expected homogeneous envelope from previous models concerning diffraction for the same 
type of experiment II16II18I . Moreover, the diffraction on the case of same branch coincidences 
is even larger than predicted by models that consider the amplitude to be constant over a slit. 




Fig. 5 Normalized two-plioton measured and simulated broad scan coincidence maps for a four- 
photon spatial state presenting a high degree of spatial correlation in each bi-photon. The upper 
row, from (a) to (c), includes the measured coincidences on the H branch (D2D4), the calcu- 
lated uncorrelated pattern coiTesponding to the product of the single counts (D2*D4), weighted 
by the repetition rate of the pump, and the SBC simulation, using Eq. U4t and Eq. )16t . The 
lower row, from (d) to (f), corresponds respectively to the coincidences between the H and V 
branches (D2D3) and to the same map corrected by subtracting the uncorrelated contribution 
(D2D3-D2*D3), as well as a simulated DBC map, using Eq. {T3j and Eq. (16). The simulation 
parameters used are, as before, a"'' = 176° and ipn''' = 170°. For a measurement time of 25 s 
per point, the maximum in map (a) corresponds to 483 counts and in map (d) to 1598 counts. 



With the purpose of understanding this, we simulated the diffraction through a single slit for 
the setup of Fig. |l(a)[ by calculating numericaUy the optical Fourier transform of the product 
of the slit transmission function and the amplitude of the bi-photon. This was done on the x+ 



and x_ directions, which are those that exhibit interesting properties for the different branch 
coincidences. We show below the optical Fourier transforms used, based once more on Eq. 
5-19of Egl: 



sinc^ (a±jc± ) oc // A(^i, ^2)^(^1, ^2)exp 

J J X[ =±A'2 



-i—[E,ixi - 
Iff 



-^2x2) 



dxidx2, (16) 



where k is the wave number of the down-converted photons, £,1 , £,2 are coordinates on the plane 
of the slits and xi,X2 are the coordinates on the detection plane. A{£,i,£,2) is the two-photon 
amplitude function on the double-slit plane and 7(^1,^2) = T {^i)T {^2) is the transmission 
function for a slit of width w, so that 7(^1,^2) = if |^i| or |^2| > w/2 and is equal to 1 
otherwise. The integrals of Eq. (fTST i were fit by the sinc^ functions and the parameters 0+ and 
a- were extracted. 

The resulting phenomenological parameters, and were used as input for the simulated 
maps presented in Fig.|5jc) and Fig.|5jf)- The simulation uses sinc^ functions onx+ andx_ with 
effective diffraction coefficients a+ and a_ as an envelope for the coincidence probabiUty maps 
of Eq. (T3[ and Eq. ( fT4] i. For the DBC case, = 1.13fl and «_ = 0.52fl, where a is the expected 
diffraction coefficient for a model that considers the amplitude constant over the length of a slit. 
For same branch coincidences, 0.29a was found by tracing over one photon of the 

pair and then proceeding with the numerical integration. We found the results to be in good 
agreement with the experimental maps. 

A similar asymmetry is reported in ll29l for a photon pairs passing through a blazed grating. 
There, also, the asymmetry is related to the degree of entanglement between the photons. 



5.2. Near-separable state 

In the near-separable case, A(^i , ^2) was indeed constant over the aperture function T{^i , ^2)- 
Thus, we found by numerical integration of the equations represented in Eq. (fTSI l that a+ ~ 
fl_ = 0.99fl. This was expected, since the photons have almost no entanglement and, therefore, 
should behave as four independent photons diffracting through the slit aperture. It can be seen 
in Fig. |6] that the envelopes for both SBC and DBC are very similar, with minor differences 
attributed in part to a small degree of entanglement present and also to image rendering artifacts. 







Fig. 6 Nomialized two-photon measured and simulated broad scan coincidence maps for a four- 
photon path state of negligible spatial congelation. The upper row, from (a) to (c), includes the 
measured coincidences on the H branch (D2D4), the calculated uncoirelated pattem con'espond- 
ing to the product of the single counts (D2*D4), weighted by the repetition rate of the pump, and 
the SBC simulation, using Eq. 1141 and Eq. (16). The lower row, from (d) to (f), con'esponds 
respectively to the coincidences between the H and V branches (D2D3) and to the same map cor- 
rected by subtracting the uncorrelated contribution (D2D3-D2*D3), as well as a simulated DBC 
map, using Eq. )13t and Eq. )16t . The simulation parameters used are, as before, a"'' = 86° and 
(p"'' = 5°. For a measurement time of 25 s per point, the maximum in map (a) con'esponds to 
75 coincidence counts and in map (d) to 1 17 coincidence counts. 



6. Discussion and conclusions 

In this work we demonstrated experimental control of quantum correlations in a four-photon 
path state comprised of two pairs of photons. Our scheme adopted a high-efficiency type-II 
parametric down-conversion source coupled to a double slit by a variable linear optical setup, 
in order to obtain spatially encoded qubits. By altering experimental parameters, we have gen- 
erated pairs with very high and very low degree of entanglement, and shown that the quantum 
correlations influence the two-photon interference pattern. That makes it possible to use a map 
of coincidence measurements for characterizing the four-photon path state. Moreover, we have 
shown that, in the case of type-II multi-pair generation, maps for coincidences between pho- 
tons with the same polarization, as well as between photons with orthogonal polarization, can 
be used to obtain a full understanding of the phenomenon without having to resort to fourfold 
coincidence detection. 

Finally, we observed that the diffraction pattem is also affected by the degree of transverse 
quantum correlation between the two photons of a pair at the double-sUt plane. We were able 
to present a semi-quantitative approach that resulted in simulated maps in good agreement to 
the experimental results. A more detailed study aiming for a closed theoretical description of 
how entanglement plays a part in distorting the diffraction patterns observed on coincidence 
measurements is under way. 

Currently, the main impediment for more efficient multi-pair experiments in quantum imag- 



ing is the fact that the area integrated by the detector is very small. The fast improvement on 
the field of CCD detectors, both in readout time and in sensitivity, opens the possibility for an 
experiment with four or more multi-pixel detectors with largely increased genuine multi-pair 
detection signal 1281 . 

Another direction of improvement is the expansion of the dimension of the logical basis from 
qubits to qudits. This can be implemented by increasing the number of slits and adjusting the 
quantum correlation control setup in an appropriate way. With a larger space, many interesting 
options become available, such as non-locality 13014321 and contextuality l33l tests. 
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